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Abstract 

In (2 + l)-dimensional QED with a Chern-Simons term, we show that spon- 
taneous magnetization occurs in the context of finite density vacua, which 
are the lowest Landau levels fully or half occupied by fermions. Charge con- 
densation is shown to appear so as to complement the fermion anti-fermion 
condensate, which breaks the flavor XJ(2N) symmetry and causes fermion 
mass generation. The solutions to the Schwinger-Dyson gap equation show 
that the fermion self-energy contributes to the induction of a finite fermion 
density and/or fermion mass. The magnetization can be supported by charge 
condensation for theories with the Chern-Simons coefficient k = Ne 2 /2-k, and 
k = Ne 2 /At:, under the Gauss law constraint. For k = Ne 2 /At:, both the 
magnetic field and the fermion mass are simultaneously generated in the half- 
filled ground state, which breaks the U(2iV) symmetry as well as the Lorentz 
symmetry. 
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I. INTRODUCTION 



Field theories in (2 + l)-dimensional space-time have been intensively studied not only 
as a laboratory for (3 + l)-dimensional field theory dynamics but also as effective theories 
at long distance in planar condensed matter physics In particular, some authors have 
investigated quantum electrodynamics in 2 + 1 dimensions (QED3) in connection with the 
effective field theories for high-T c super-conductivity [fj].[] In this theory, there can be a 
topological, i.e. metric independent, gauge action, known as the Chern-Simons (CS) term, 
which connects a magnetic field B with an electric charge density e(ip^ip) for the fermion 
field ip. Because of this peculiar property, the CS term is essential in the field theoreti- 
cal understanding of the fractional quantum Hall effect, which is based on the composite 
excitation of an electron and magnetic fluxes 

A theory whose gauge field action includes both the CS term and the Maxwell term 
was proposed by Deser, Jackiw, and Templeton [[7] as an attempt to improve the infrared 
photon behavior in QED 3 . In this theory (Chern-Simons QED 3 ) the coefficient k for the CS 
term gives the photon a gauge invariant mass which explicitly violates parity. Hosotani || 
showed that spontaneous magnetization occurs in Chern-Simons QED 3 , because the gauge 
invariant photon mass k is completely screened out by parity- violating vacuum polarization 
effects. This also means that Lorentz symmetry is broken through the induced magnetic 
field and the resultant massless photon plays the role of the Nambu-Goldstone (NG) mode. 
In this theory, the Gauss law, 

«S = -e(VV), (1-1) 

follows from the equation of motion. Thus the magnetized vacuum corresponds to the ground 
state with finite fermion density, ("0 V) 7^ 0, that is, an occupied lowest Landau level. 

In Ref. ||, however, the finite density vacuum configuration was set by hand as a state 
with a definite filling fraction v that represents the degrees of occupation of fermions in the 



lowest Landau level. The Gauss law (|1 . 1| ) tells us that spontaneous magnetization means a 
phase transition occurs between different vacua assigned to specific values of v. However, 
without any parameter connecting differently occupied vacua, it was not clear whether the 
condensate (if) >i[>), which supports the magnetic field, is spontaneously induced or finely 
tuned. 

In order to clarify the above situation, an external source term fiifj^ip was introduced 
in Ref. to control the finite densities in the vacuum with a continuous parameter /x. 
From the viewpoint of statistical mechanics, /1 plays the role of the chemical potential. The 
parameter \x has to exceed the fermion mass m, or the energy of lowest Landau level, in order 
to induce a finite fermion density. We notice, however, from the field theoretical viewpoint, 
the chemical potential term /j ip^ip also plays the role of a term that explicitly breaks Lorentz 



x The application of gauge field theories to high-T c super-conductivity was originally proposed 
and studied in Ref. || in the context of the resonating valence bond state. The non-Fermi liquid 
behavior in the normal phase of the high-T c system has been studied with the help of effective 
gauge theories and their renormalization group analysis. See Ref. [|| for details. 
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symmetry. So long as the fermion mass m remains constant, we cannot take the Lorentz 
symmetric limit fi — > while retaining non-zero finite density in the vacuum. Alternatively, 
however, if we take the symmetric limit such that m < |//| — > 0, the finite density vacuum 
might be spontaneously realized through the condensate (ip'i/)). In other words, the finite 
density vacua are spontaneously realized if and only if the fermion bare mass becomes zero 

Recently, in Refs. [0,0] it was shown that a strong magnetic field enhances the conden- 
sate (V'V')j or the dynamical generation of fermion mass, in (2 + l)-dimensional four-fermion 
interaction models. Refs. P"U|,|T2"| explained this effect as the dimensional reduction of the 



phase space for charged particles coupled to a strong magnetic field.f] In a strong magnetic 
field, the wave function for charged fermions is localized within a region whose size is given 
by the magnetic length: 1 = 1/ \J\eB\. The fermion in 2 + 1 dimensions, therefore, behaves 
like that in + 1 dimension, while the photon field is charge neutral and propagates freely 
in 2 + 1 dimensional space-time even in the presence of the background magnetic field. The 
condensate (ipip) is therefore easily formed much like in Bardeen-Cooper-Schrieffer (BCS) 



theory |I|[TJ. 

The Lagrangian of QED 3 with iV flavor massless four-component fermions has U(2iV) 
symmetry. When N is smaller than its critical value 3 < iV c < 5, the flavor U(2iV) symmetry 
is spontaneously broken by the condensate (ipip) and the fermion acquires a dynamically 
generated mass |15|JT6|] . By means of the Schwinger-Dyson (SD) gap equation, Shpagin 



|T7| showed that a fermion mass is dynamically generated and the U(2iV) symmetry is 
spontaneously broken irrespective of N in QED 3 with an external magnetic field. 

This result means that the magnetic field catalyzes the fermion mass generation and 
seems contradictory to the results of Refs. [§,|9] at first sight, since the magnetic field is 
spontaneously induced only for massless fermions in Refs. |8|||. If both results of Refs. 
HU and Ref. [0 are true, the chemical potential fiip^ip should affect the condensate {ipip) 
in the presence of a background magnetic field. It has obviously become an interesting 
question in Chern-Simons QED3 to study the dynamical generation of both condensates, 
(iji/j) and V)> or the fermion mass and the magnetic field, as a specific (2 + l)-dimensional 
phenomenon. 

In this paper, we will directly show that the spontaneous magnetization, or Lorentz 
symmetry breaking, occurs as the realization of the finite density vacuum in Chern-Simons 
QED 3 . We will also examine the possibility that Lorentz symmetry is broken in a theory 
where both massive fermions and a magnetic field are spontaneously generated [18 . 



The paper is laid out as follows. In Section [TT] we derive the effective action in the large N 
limit that contains the screening of the iV flavors of fermions. We add to the Lagrangian both 
the fermion bare mass term mipip as an explicit breaking source term for the flavor U(2iV) 
symmetry, and the chemical potential term fiip^ip for Lorentz symmetry. In Section [1 1 1 
we directly calculate the condensates, (^0) an d as the large N quantum correction 

of the fermion field in the presence of the background magnetic field. The effect of the 
explicit breaking parameters and the magnetic field B are taken into account in the fermion 
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The enhancement of fermion mass generation by a magnetic field was found also in lattice QED3 
[13] and was applied to the effective field theories in planar condensed matter systems [ |l4|| . 
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propagator by using the proper time method |L9| . Thus we can see the dependence of m and 
/i upon the condensates, (ipip) and (ip^ip), explicitly. It is found from the Gauss law (|1 . 1| ) that 
the consistency condition for nonzero magnetic field places a restriction on the CS coefficient 
k in the symmetric limit fi — > 0, m — > 0. In Section [TV] we calculate the vacuum polarization 
tensor in the large N limit in a gauge invariant manner. We regularize the parity-violating 
part, or the induced CS term, in a manner consistent with the charge condensation. The 
improved photon propagator including the effect of vacuum polarization is also derived. In 
Section [V| we construct and solve the Schwinger- Dyson (SD) gap equation for the fermion 
self-energy that contains the dynamically induced fermion mass and chemical potential 
fid- We investigate the solutions to the SD equation in the symmetric limit fi — > 0, m — > 0, to 
clarify whether or not the fermion mass and/or the magnetic field are dynamically generated. 
Section [VT] is devoted to the study of vacuum stability. We calculate the Cornwall- Jackiw- 
Tomboulis (CJT) potential |2(J for the composite operators, ipif) and if)' if), and the shift of 
zero-point energies for the fermions and the photon. In Section |V11| we conclude our paper 
with a discussion of the relation between the results in Refs. |8]||, UTTH , and ours. 

II. LARGE N EFFECTIVE THEORY IN CHERN-SIMONS QED 3 

In this section we construct the large N effective action in Chern-Simons QED 3 , which 
includes long range screening of charged fermions. We use the metric diag(g M1/ ) = (—1, 1, 1) 
through out the paper. 



A. The model and the symmetries 

The Lagrangian density of Chern-Simons QED3 with massless fermion is 

£ = ^ [id, + eA,} if) - -^F" - ^ p A»d»A p - ^{dAf, (2.1) 

where the third term is the Chern-Simons term which is topological (or metric independent) 
and the last term is a covariant gauge fixing term. From the Lagrangian (|2.1| ) we obtain the 
equation of motion 

dvJW - ^ vp T vp = -eif)^if) (2.2) 
whose vacuum expectation value leads to the Gauss law constraint kB = —e(if)'if)) [cf. Eq. 



(|1.1|)1, under a given magnetic field. 

The smallest spinor representation of Lorentz group in 2 + 1 dimensions is provided by 
two-component spinors. A corresponding 2x2 matrix representation of the Clifford algebra 
is expressed by the Pauli matrices 

7 = 03, 7 1 = ia x , 7 2 = ia 2 , (2.3) 

which obey the anti-commutation relation: {7^,7^} = —2g fiu . We have no other 2x2 
matrix which anti-commutes with all 7^. We, therefore, have no chiral symmetry that 
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would be broken by a mass term rmpip in two-component representation, but the discrete 
parity symmetry is broken by the same mass term. To construct a continuous group which 
resembles the chiral symmetry in 3 + 1 dimensions, it is necessary to enlarge the spinor 
representation to be the one with four-component spinors. Now we have a Clifford algebra 
given by three 4x4 matrices, namely, 

We can then construct another two 4x4 matrices 

7^(° fU = *f° i), (2.5) 





which anti-commute with 7 , 7 1 , and r y 2 . Using 7 s and 7 5 , we can define a U(2) group with 
generators 

t = I, n = 7 5 , r 2 = -ry 3 , and r 3 = 7 3 7 5 . (2.6) 

We see that the Lagrangian density for a massless four-component fermion is invariant under 
this global U(2) symmetry. The appearance of mass term mipip breaks the U(2) symmetry 
to a subgroup U(1)®U(1) whose generators are r and r 3 . For the system of N flavor four- 
component fermions, global symmetry is extended to U(2iV) symmetry which is nothing but 
the direct product of the U(2) symmetry and the flavor U(iV) symmetry. Once again, the 
mass term rmpip breaks this U(2iV) symmetry to its subgroup U(iV) ® U(iV). 

In 2 + 1 dimensions, the inversion of two space axes is equivalent to a ir rotation on the 
two-dimensional space so that parity transformation should be defined as an inversion of one 
spatial axis: P, (x,y) — > (x,y)p = (—x,y). If we define the parity transformation of field 
operators, it should leave the Lagrangian of massless QED 3 remain invariant. A specific 
construction is to transform a two-component spinor and the gauge field as follows; 

P^(t,x)p- 1 =<7iV>(t,X P ), 

PA^(t,x.)p- 1 = (-l)^^(i,xp). (2.7) 

A significant feature of the two-component Dirac spinor is the fact that mass term rmpip 
changes its sign under the parity transformation P. As to the gauge field the CS term is 
odd under the above parity operation. It has been known that a two-component fermion 
induces a CS term with the same signature as the fermion mass through the one-loop vacuum 
polarization [[7,21 . 



A convenient method to construct the parity even mass term for fermions in 2 + 1 
dimensions is to use a four- component fermion in terms of two two-component fermions ipi 
and ip2 '■ 



*=^J- (2 ' 8) 

Parity transformation for the four- component spinor ip is defined as 
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Pip(t, x)P _1 = -i7 3 7 V(^, x P ). (2.9) 

It includes the exchange of upper and lower two-component spinors, ipi — > ip2, ^2 — * in 
addition to the P operation for each two-component spinor. Therefore it is easy to write 
the parity-conserving mass term: 

mipijj = m^j\a 3 -ipi — rmpla 3 tp 2 , (2-10) 

that breaks the U(2iV) symmetry to the U(iV) <8> U(N), as well as the parity- violating mass 
term: 

that is a singlet of the U(2iV) group. 

According to the theorem found by Vafa and Witten []22 1 , it is energetically favorable that 
the U(2iV) symmetry is broken to the vector-like global symmetry U(iV)®U(iV) in the vector- 
like gauge theories such like QED3. On the other hand, the theory of our interest contains the 
topological Chern-Simons term as well as the background magnetic field. Since the theory 
ceases to be vector-like in general, the stability of U(iV) <£)U(iV) is not necessarily supported. 
In this paper we constrain ourselves to the case that the U(iV)®U(iV) symmetry is unbroken, 
therefore the fermion mass is dynamically generated only as the parity-conserving mass term 
and the induced CS term due to the fermion mass is completely canceled. 



B. Effective action in the large N limit 

It is known that QED 3 is a super-renormalizable theory and its /5-function for a dimen- 
sionless coupling constant, which is defined as the photon propagation in the large N limit, 



shows a nontrivial infrared fixed point |]15 |. On the fixed point the gauge coupling e scales 
like e ~ 1/ y/N, which means that low energy effective theory is dominated by large N con- 
tributions. In other words we should not neglect the screening due to vacuum polarization 
at long distance. In the following, we construct the effective action in the large N limit, 
which rules infrared dynamics of Chern-Simons QED3. 



According to Ref. |TJ|, we introduce a dimensionful coupling a = Ne 2 /Ait and keep a 
finite when N is taken to infinity so that radiative corrections can be added to the effective 
action successively in 1/N expansion. What we want to ask is whether or not a magnetic 
field is dynamically induced. On the ground state with the induced magnetic field, proper 
excitations are the fermions coupled with the magnetic field and the gauge field fluctuation 
around the magnetic field. Taking into account the above expected physics, we divide the 
gauge field A^ into a background field and a fluctuating field; A^ = y/N (A^ + A^) with 
\/WA e * t (x) = Bx^fjix. Notice that we have rescaled the gauge field by a factor yN, which 
enables us to treat loop expansion as 1/N expansion. Subsequently, the magnetic field B 
must scale as y/N in order that eB remains finite in the large iV limit. 

The Lagrangian density ( |2.1|) is divided into the gauge part and the fermion part; £ = 
£>g + £>f, 

£g = - ^e^A^A" - ^{dAf - ^ - VNkBA , 

C F = tfj lij^DJA™ 1 + A] -m + ^7°! V, (2.12) 
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where -D/JA] := — ieyNA^ is a covariant derivative. The propagator for the fermion 
coupled to A° xt is defined by 



S(x,y) 



x 



1 



, r^Z^A 6 **] -m + /i7° 
It is written in the proper time method |0| such as 

'ie 



S(x, y) = exp i^-VN (x - yYAf{x + y) ) S(x - y). 
Fourier transformation of S(x — y) is given by 

S{k) = i 



(2.13) 



(2.14) 



roc 




/ ds exp 


—is ^ 


h 





2 _ k 2 + tan(e£s) k2 



x { [l + 7 V tan( eJ Bs)] (m + 7%) - (7% + 7%) sec 2 (eBs)} , (2.15) 

where k e := k° + fi + iesgn(k°) modifies the ie prescription to be consistent with the shift 
of Hamiltonian by \i [24|. In Eq. Q2.15| ) the integration over s is the only formal expression 
and, rigorously, we have to choose the sign of s according to the sign convention of k° + fi 
so that the integral over s converges. The photon propagator from Lq is 



1 



p 2 + K 



p^p , . p P t 

5- + IK 



livp 



P 



,2^2- 



(P 



(2.16) 



which tells us that the CS coefficient k behaves like a pole mass of the gauge field. 

The large iV effective gauge action is obtained by integrating out the fermion field |25 
that is, 



T[A] = -An / [#][#] exp 



Jd 3 x £ 



d 3 x C G - iN Tr Ln -i^D^A 6 ** + A] + m - /i 7 c 



(2.17) 



To expand Tr Ln term with respect to the gauge field A^, we will extract from it charge 
condensation, vacuum polarization, and nonlocal self-interactions of the gauge field. Specif- 
ically, it is 



T[A] 



-iN Tr Ln S" 1 + N J d 3 x 
+ N [d 3 x 



B 



K B 1 

— - -^=A°(x) - -A»(x)A~l{-id x )A u (x) 



N 



A" x) 



d 3 x A"{x)U l)V (-id x )A v (x) 



+ f[Aj, 



(2-18) 



where (ip'ip) is the charge condensation induced by a magnetic field and obeys the Gauss 
law. T[A] denotes the nonlocal vertices of the gauge field. Vacuum polarization tensor 
is given by 



n™(p) = -Ne 2 



d 3 k 
i(2vr) : 



tr 



lnS{k)~f u S(k-p) 



(2.19) 
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which improves the kinetic term of the gauge field. 

By introducing the improved photon propagator: D^ u (p) := [A~J(p) — n /u/ (p)] -1 , the 
above effective action in a complicated form is simplified as 

F[A) = -iNTiLnS- 1 -il^-- y J^xA^x)D^(-id t )A v (x)+f[A], (2.20) 



where Q is three-dimensional space-time volume. The linear terms for in Eq. ( 2.18|) have 



been canceled out by means of the Gauss law. Later in Section [V] we will study dynamical 
symmetry breaking through the analysis of SD gap equation for the fermion self-energy. It 
can be derived from the effective action ( |2.20| ). In Section |Iy| vacuum polarization tensor 
H /Jil/ (p) is regularized in a gauge invariant manner, such as 

iWp) = (p^v - v 2 g^) n e (p) + (p>^ - p\g^ v ) n ± (p) - i P % up n G (p), (2.21) 

where p^_ = (O,^ 1 ,^ 2 ) and diag(g^) = (0, 1, 1). The parity-violating part U. (p) provides a 
quantum correction to the gauge invariant photon mass k in the propagator D^ u . It will be 



explained in Section |T| that the Gauss law forces the improved photon mass, k — II o (0), to 
vanish in the presence of a magnetic field. Therefore our effective theory in Eq. ( |2.20| ) has 
the massless photon, whose propagator is given by D^ u , as well as the low energy effective 
theory in QED 3 with a magnetic field. 



III. CONDENSATES AND CONFIGURATIONS OF VACUUM 

In this section we explicitly calculate two condensates, (ipip) and ("0 V)' f° r fermions 
coupled with a magnetic field. We can find the possible configurations of vacuum (the lowest 
Landau level) by taking both condensates to the symmetric limit J± := fi ± m — > 0. The 



condensates are calculated from the fermion propagator ( |2.15| ) which contains the explicit 
symmetry breaking parameters, fi and m. We shall find various patterns of symmetry 
breaking according to different approaches to the symmetric limit J± — > 0. It is also clarified 
by applying the result of calculation to the Gauss law that possible values of the CS coefficient 
k are restricted to \k\ =0, a, and 2a. 

The propagator for iV flavor fermions ipi (i = 1, ...,N) coupled with a magnetic field is 
defined by 

S(x,y)8 ij :=z(0|T^(x)^-(y)|0), (3.1) 

where we assume the U(iV) ® U(iV) symmetry and extract the flavor indices as a diagonal 
metric 5^. S(x,y) in Eq. ( |3.1| ) is equivalent to the definition in Eq. ( ^.13[ ). We introduce 
the functions, J and Jq, which are defined by the following equations: 

J(m,fJL) := - lim tr (0|T^(x)^(y)|0) 
y 

= iNtr S{x,x), (3.2) 

J (m,n) := - lim tr 7 °(0|T^(x)^(y)|0) 
y >x 

= iiVtr j°S(x,x). (3.3) 
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The condensates, (ipip) an d {ip , are identified as J and j7o in the symmetric limit J± — ► 
such as 

(^) = lim J( m ,n), (V'V) = I™ J (m,^). (3.4) 

If the condensate (ipip) (('0'V'}) remains nonzero, the vacuum is realized as a non-singlet 
state of U(2iV) symmetry (Lorentz symmetry). Precisely, the symmetry is spontaneously 
broken by the condensate though the U(2iV) symmetry (Lorentz symmetry) is restored at 
the tree-level Lagrangian in the symmetric limit J± — > ||. 



A. Calculations of condensates 

Since the condensates will be computed in terms of the fermion propagator ( |2.15| ), a 
comment about the meaning of Wick rotation, k° = ik, in the presence of [i should be noted. 
The improved ie prescription in Eq. ( |2.15| ) always forces poles of the fermion propagator on 
a complex k° plane to be located outside the contour which connects the real axis to the 
imaginary axis. Therefore the Euclidean propagator becomes 



S(k E ) 



ds exp 



-s \m 2 + (k - %[if + tanlir fc2 

0/ 



| 1 — ry^ 2 sgn(e.B) tanhr m + ry°(A; — z/i) — ^y±k± sechVj 



(3.5) 



by using r := \eB \ s. We will carry out analytic continuation k — > k — ifi after integrating 
over a proper time s. Hereafter the shifted momentum k — ifi in the proper time integration 
necessarily implies the above agreement even if we do not indicate it explicitly. 

Now we calculate the condensates J and j7o- They are written in terms of the Euclidean 
propagator S{k,E) as follows; 



J{m,fi) = -N 



d 3 k f 



tr S{k L 



(2tt) 3 

N 



Jo(m,/j) 



2tt 2 
-N 



\eB \ ml dk ds e 

J-oo Jl/A 2 

<H '" tr ^S(k E ) 



s[(fc— i^i) 2 +m 2 



coth r, 



{2*y 

eB i 



2tt 2 



dk (k — i/j,) / ds e 
Jo 



coth r, 



(3.6) 



(3.7) 



where the Gaussian integrals for k have been done and a cut-off scale A has been introduced 
to regularize the ultraviolet divergence in J7"(m, fi). To extract finite density effects, we first 
decompose the integrands into Landau levels by using an identity: 



cothr 



n=l 



-2nr 



(3.8) 



Then we integrate out the proper time s and get 



S 



N 

J(m, /i) = J(m, 0) - — — \eB\ m 

Z7T 



/ 1 °° / 1 

Jc„ ^ z 2 + m 2 + JCp ^ z 2 + £ 2 



J Q (m,ix) 



N 



2tt 2 



\eB\ 



dz 



iz 



Cu z 1 + m' 



+ 2E 



dz 



'IZ 



i J c» z2 + Q. 



(3.9) 
(3.10) 



where £ n := ^m 2 + 2n\eB\ is the n-th Landau level. The above integrations are carried 
out along a contour on a complex z-plane depicted in Fig. [I]. The contribution at zero 
density appears only in the condensate (ipip) through J7"(m, 0) which is given by 



J(m,0) 



N 



2tt 3 / 2 



\eB \ in 



00 ds 



l/A 2 s 



1/2 



coth r. 



(3.11) 



Using contour integrals shown in Appendix 0, we obtain the following results; 



J{m,fi) 



^"777 Am \eB\ 

7T 3 / 2 27T ' ' 



/- , ( I (ml) 2 
1 + V2ml( 



iV 





%|-») + 2EfW-^) 

n=l 



Jo(m,fi) = —1651 sgn(/i) 

Z7T 



n=l 



(3.12) 
(3.13) 



where we have summed up the Landau levels with the generalized Riemann zeta function 
and the value of step function 9{x) at x = is chosen to be 1/2 at zero temperature limit 
of Fermi-Dirac distribution. 

We notice that Jo reduces to zero as we take the limit, B — > (I — > oo), in Eq. ( [3.131) . 



The charge condensation J§ is known to be unaffected by higher-order corrections [26 



Therefore the above one-loop result is exact as we shall make detailed comment later on in 
this section. 



B. Configurations of vacuum 

It should be noticed that we can read a filling factor v as a function of \x from Eq. (|3.13|) 



26|j27[| . Since the wave function for fermions in the lowest Landau level spreads out within 



a size of the magnetic length 1—1/ J\eB\, the number of states for an electron is estimated 
as l/27r/ 2 per unit area. We also know that /i > (// < 0) induces a finite density of fermion 
(anti-fermion) whose spin indicates the up- (down-) state in the lowest Landau level. 
is rewritten as 

(^) = £p{^+-^-), (3-14) 

where u + (z/_) denotes the filling factor of fermion (anti-fermion) whose spin indicates the 
up- (down-) state. Comparing Eq. (|3.13| ) with Eq. (|3.14|) , we read 



v± = e(±n)B(\n\-m). (3.15) 
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Three different approaches to the symmetric limit are possible on a m) plane, namely, 
\fi\ < m — > 0, m < — > 0, and = m — > 0, in which both condensates, (^V) an d (^ip), 
exhibit the different patterns of symmetry breaking (see Fig. |2p 15). 

For < m — > 0, we approach to the origin of the (|/i|,m) plane along a line inside the 
upper area: < m. The condensates can be read as 

W> = -^ (^V) = 0, (3.16) 

which shows spontaneous breaking of the U(2iV) symmetry on an empty vacuum, that 
stands for the lowest Landau level with v± = 0. If we believe the Gauss law (|1 . 1|) naively, 
the generation of magnetic field is available only when charge condensation occurs. Since 
we have no charge condensation seen in Eq. ( |3.16| ), we expect no spontaneous magnetization 
and, subsequently, no spontaneous breaking of the Lorentz symmetry. However, we want to 
postpone our conclusion until we study the vacuum energy shift due to the magnetic field 



in Section VI 



For m < \fj,\ — > 0, we choose a line inside the lower area, m < |//|, to approach the 
symmetric limit. We get an opposite pattern of condensates to the former case, namely, 

iV 

(W) = 0, {^) = —^aga(ji), (3.17) 

which mean the vacuum, or the lowest Landau level, has been filled up as it is read from 
(z/ + ,i/_) = (1,0), or (0,1), and the U(2iV) symmetry has been restored. This time we 
have obtained a finite density vacuum so that the magnetic field is supported by nonzero 
charge condensation. Therefore the magnetic field can be generated dynamically. If the 
gap equation for fermion self-energy favors a finite density solution, Lorentz symmetry is 
spontaneously broken through the induced magnetic field. 

For = m — > 0, the symmetric limit is taken along the boundary line of |/i| = m on 
the (|/i|,m) plane. Both condensates remain finite and are given by 

<^> = -^, (V>V> = sgnGu), (3.18) 

which mean that the U(2iV) symmetry is broken spontaneously in the half-filled vacuum 
with (z/ + ,z/_) = (1/2,0), or (0,1/2). As to the Lorentz symmetry, it might also be broken 
spontaneously through the induced finite density which supports the magnetic field through 
the Gauss law. 

We notice that the two condensates always appear in such a complementary way as 
to keep the combination {ipip) — \ (ip^ip)\ as a nonzero constant: —N/2nl 2 . Moreover the 
above result depends only on the kinematics of fermions in the lowest Landau level. Even 
if we have no attractive interaction of photon exchange, the magnetic field alone can cause 
the nonzero condensate: (ipip) — \(i[>>i{))\. This is a specific phenomenon realized only in 
(2 + l)-dimensions, and, for (^ip), the same discussion was already made in Ref. [|Kj. 



3 The restoration of U(2iV) symmetry in a given finite density was also pointed out in Ref. [1C] 
in (2 + l)-dimensional Nambu-Jona-Lasinio model with a magnetic field. 
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C. Consistency condition 



In Chern-Simons QED 3 , charge condensation {^i/j) is necessarily connected with a mag- 
netic field through the Gauss law ( |1.1| ), while Eqs. ( p. 16| - p7T^ ) tell us that the magnetic field 
keeps (ip'ip) as quantized values in a unit of l/27r/ 2 even in the symmetric limit. Therefore, 
the Gauss law places a restriction on the CS coefficient k in order to retain the magnetic 
field in the symmetric limit. Substituting Jo(m, /x), given by Eq. (|3Tl3|) , for (ipty) in the 
Gauss law (|1.1|), we obtain in the case of < Si the equation 



B 



n + 2a sga(iieB) 6(\n\ - m) = 0, (3.19) 



which is regarded as the self-consistency condition for a given magnetic field B. In order to 
keep the nonzero magnetic field in our system, k should be restricted to 



K 



' for < m — > 

-asgn(fieB) for = m — > . (3.20) 
-2asgn(/ie5) for m < — > 



If we start from the theories with \k\ = 0, a, and 2a, which allow the penetration of 
nonzero magnetic field, then each symmetric limit (|/i|,m) — > (0,0) is uniquely determined 
as |//| < m — > 0, = m — > 0, or m < — > 0, respectively. 



The above constraint for k is also interpreted as follows. It will be shown in Section [TV 
that, except for = £ n (n = 1, 2, . . . ), the charge condensation Joijn, /i) is related to the 
parity- violating part of vacuum polarization II (p) through the relation |2 



-ej {m,n) = BU o (0). (3.21) 

The constraint ( |3.19| ) is interpreted as B[k — n o (0)] = which means that a magnetic 
field applied to the system can penetrate it if and only if the effective photon mass, 
k cS := k — II o (0), becomes zero P J2S|J29] ]. Otherwise, the system with k ^ n o (0) ex- 
cludes a magnetic field irrespective of its origin. In this case we should set magnetic field 
zero and investigate dynamical generation of the parity- violating fermion mass as well as 
U(2iV) symmetry breaking [|30| , [31|1 .fl In the theories with k = Il o (0), we cannot exclude the 
possibility of spontaneous magnetization since we achieves the symmetric limit which allows 
the penetration of a given magnetic field in each theory [cf. Eq. ( |3.20j )1. We notice that 
the above statements are not affected by any higher-order correction by virtue of the non- 
renormalization theorem of the induced CS coefficient IT o (0) (Coleman-Hill theorem p2"| ) 
even with finite density and magnetic field |26],|33 . 



In Section [V] we investigate dynamical generation of fermion mass and magnetic field 
for theories with \k\ = 0, a, and 2a. We shall confirm that the fermion mass and the 
chemical potential /i^ are dynamically generated in the symmetric limit (|/i|,m) — > (0,0). 
In order to retain the consistency with a nonzero magnetic field, and Hd should satisfy 
the same patterns as the explicit breaking parameters m and /i do in Eq. (|3.20|) for each k. 



It was pointed out in Ref. [31] that spontaneous breaking of the U(2iV) symmetry turns out to 
be a first-order phase transition in the theories with k ^ n o (0). 
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IV. VACUUM POLARIZATION 



This section is devoted to the gauge invariant derivation of the vacuum polarization 
tensor U^ v defined by Eq. ( j2.19|) , that is 

— ^ tr [ ltl S(k) lu S(k - p)] . (4.1) 

It should be regularized so as to be gauge invariant such as 

IV(p) = (PuPv - P 2 9^) n e (p) + (p^pi - p\g^) n ± (p) - ip p e livp U {p). (4.2) 



As mentioned in Section |TJ, effects of an explicit breaking parameter /i, or equivalently a 
chemical potential, are included through an analytic continuation, k — > k — ip, of the third 
component of an Euclidean loop momentum. As far as we derive parity-conserving parts of 
vacuum polarization, they are obtained through the above analytic continuation from zero 
density results. In the following we divide the vacuum polarization tensor into zero density 
parts and finite density parts, that is, 



U e (p) = W e (p) + U{(p) 

rf 



n ± (p) = ni( P ) + nl(p), (4.3) 
n (p) = ul(p), 

and derive them separately. In our setting of N flavor four-component fermions the parity- 
violating part n o (p) appears only as a finite density effect f\ and is uniquely determined 
through a gauge invariant constraint between II o (0) and charge condensation. 



A. Gauge invariant regularization 

The regularization of U pu is worked out in the same gauge invariant manner as in the 
case without magnetic field. First we set p = and mix two proper times of S(k) and 
S(k — p) by means of a parameter integration. Then we carry out the Gaussian integrals for 
spatial components of the loop momentum: k. We find that parity-conserving parts appear 
only as even terms for k° and include ultraviolet divergences, while a parity-violating part 
appears only as a linear term for k° and is finite. 

As to the parity-conserving parts we see that the ultraviolet divergences appear as terms 
with the lower power of proper time s, namely, s -3 / 2 in 2 + 1 dimensions. Integrating out the 
proper time s by part, we can subtract the ultraviolet divergences as boundary values of the 
integration over s. Then we find that remaining finite contributions satisfy the transverse 
form in Eq. (|4.2|). Zero density parts of IF and IF are determined as |23| 



Since in zero density a parity-violating part necessarily appears as an integration of an odd 
function of k° which is the zeroth component of the loop momentum, it vanishes after an integration 
over k°. In our setting of fermions the parity violation due to the fermion bare mass m is completely 
canceled out. 
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Kip) 
ni(p) 



a 



2^J- 
a 



f 1 d f°° ^ S - s [ m2 +^(p)] T ( cos ^ TV ~ v cothrsinhrf) 
7-i Jo s 1 / 2 sinhr 

J -i Jo s^ 6 sinh 3 r APh 



(4.4) 



where we have introduced a rescaled variable r := s/~ 2 and the function <p(p) is defined by 



<f>(p) :-- 



1 — v 2 2 cosh r — cosh tv 2 



4 



-P\\ + 



2t sinh r 



P_l- 



(4.5) 



In order to derive the finite density contributions to the parity-conserving parts we only 
have to insert Gaussian integral: 



dk e 



1 



7T J-oo 



(4.6) 



into the expressions of Ilg and H z ± in Eq. ( f4.4| ) and to perform analytic continuation k — > 
— ifi. 

As to the parity- violating part it is given as a combination of SO (2) invariant tensors 
under a given magnetic field, that is, 



n^(p) 



a 

2^ 

x 



dk(k-ifi) f 1 dv [°° dse- s l Ck -^ 2+m2+ ^ 



2r(cosh r cosh tv — 1) p / y j_ || j_ \ 

P± (5^0 e i/p — Quo e n P ) 



sinh t 



2tv sinh rt> n 
+ — P°e 



sgn(e J B), 



sinh t 

which vanishes for \x = after an integration over k. Even if we use an identity 







V Pt nvp — P± [Q'/lo e U p 9uo t 



(4.7) 



(4.8) 



in order to extract the gauge invariant tensor p p e ilvp , there remains an ambiguity in n o (p). 
However, in the presence of the magnetic field, there exists a gauge invariant constraint 
between n o (0) and charge condensation eJio(m, fj,) f26|, namely, 



dJo(^m) 

By virtue of Eq. (|4.9|), II (p) is uniquely determined from the result of Eq. ( ft.lcp . 
Through the above procedures, the functions ITg, IT^, and n„ are determined as 



(4.9) 



n f e (p) 



ni(p) 



a 
2^ 



dk dv ds <e 



-s[(k-ip) 2 +m 2 +4>(p)] _ -s[k 2 +m 2 +<f>(p)] 



r (cosh tv — v coth r sinh tv) 



a 
2^ 



sinh r 

1 PO 

dk I dv I ds \ e 



i[(k-ip) 2 +m 2 +<t>{p)\ _ e -s[k 2 +m 2 +(f>(p)] 
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2r(cosh r — cosh tv) _ f , . 

x -T3 K(P)> 

sinh r 



(4.10) 



dkiik - ifi) J 1 dv J™ ds e -4(fc-^) 2 W+0( P )] 



a 
2^ 

2r(cosh r cosh tv — 1) . _,. 

x —2 sgn(eS). 

smn r 



As mentioned in Section PTT] , one cannot do the analytic continuation k k — ifi before the 
proper time s is integrated out, and then the order of integrations should be kept in Eq. 

When the momentum p vanishes, it is possible to integrate out the parameter v and the 
proper time s before the integration over k. Therefore we can look into how the analytic 
continuation works. The integrands in Eq. ( 4.10|) are decomposed into the Landau levels by 
means of identities: 



r (cosh tv — v coth r sinh tv) 
sinh r 
2r(coshr — coshrf) 
sinh 3 t 
2r(cosh r cosh tv — 1) 
sinh 2 r 



r 



oo 

5^[l + (2n + l)v 



n=0 
oo 



e - { 2n + l + v)r + , _s 



At [n 2 e' 2nT -n(n + l )< 



n=l 
oo 

-2r 



-(2n+l+u)T 
2nr /o„ i i N „-(2n+l+v)T 



(4.11) 



^ [2ne~ 2nr - (2n+ l)e 



n=0 



+ (v 



Inserting the above identities into Eq. ( |4.4j ) at zero momentum (p = 0), we acquire the 
following results for zero density parts after we integrate out the parameter v and the proper 
time s; 



II*(0) = -2al 
ni(0) = 2al 



r ( 1 (ml)' 
ml + 3V2( T2' 2 



V2 



1 (ml)' 



+ 1 



mZ H — — ( 



1 (mZ)' 



V2 



2' 2 



+ i +^-K( - 



3 (ml)' 



+ 1 



(4.12) 



where we have summed up contributions from the higher Landau levels by means of the gen- 
eralized Riemann zeta function. Substituting Eq. (^4.11|) into Eq. fl4.10|) at zero momentum 
(p = 0), we also obtain the following expressions for the finite density parts; 



n f e (0) = -^j£ dzl 2 \nl 2 (z 2 + m 2 ) + 2j^ £ ^Z 2 ln, 

a 

71 



ni(o) = ^£ 



l 2 (z 2 + S 2 ) + 2n 



dz 



n=l 



-n I dz — — + — / dz 

Jc, z* + ei p Tc (z 2 + ^) J 



1^(0) = — 8gn(e5){ / + fdz 

7T { Jc^ z 2 + m 2 ^ [ Jc„ 



lc„ z 2 + £, 



K(0); 

2n 
' ~ I 2 



c, z 2 + S 2 _ 

(4.13) 



dz 



iz 



[z 2 + S 2 ) 



When = £ m poles or branch points on a complex z-plane are located on the contour 
and force us to define the integrations as their Cauchy's principal values. Subsequently, we 
obtain (see also Appendix |A|) 
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{oo 
l(\fi\ -m)e(\n\ -m) +2^ 



n f ± (0) = 4al^T 



n=l 



2n 



rr 



lE n (lE n 



n=l 



l(\fl>\ - £n 



n 



n 



l(lS n 



6(\»\-E n ) -n f e (o), 



(4.14) 



II f o (0) = -2asgn(/ie J B) 9(\fi\ - m) + 2 ]T 



n=l 



<5(|£i| - E n ) 



Now we prove Eq. ( |4.9| ) that is the gauge invariant constraint between the charge conden- 
sation and the parity- violating part of vacuum polarization and is essential to determine the 
n o (p) uniquely. We see the following relation is obtained from the definition of the fermion 
propagator [cf. Eq. (|2.13|) ]; 







OB 



—S{x, y) = / d 3 z S(x, z)VNe ^ — ^-^S(z, y) 



dA^Hz 



OB 

e \ d 3 z Z2 S(x, z)^ 1 S{z 1 y), 



(4.15) 



where we have used the explicit form of the gauge potential v4^ xt (z) in the second line. 
Applying the above equation to the definition of charge condensation in Eq. (|3.3|), we obtain 



d 

e 'Qg<Jo{ m ' A 4 ) = iNe 2 1 d 3 z z 2 tr j°S(x, z^Sfe, x) 



d 3 p 5 3 (p) (x 2 - U 01 (p) = -n c (o). 



(4.16) 



We see that in Eq. ( [4. 14j ) delta function singularities appear only in the higher Landau 
levels, E n [n > 1), which have gap energies owing to the magnetic field B. In order to take 
the symmetric limit J± — > 0, it is enough to constrain the region of parameters, J±, to be 
on the lowest Landau level. For |/i| < E\ we notice that Eq. fl4.9|) reduces to 



eJoK/i) = sn o (0), 



(4.17) 



which has been used in Section [TTI| to derive the consistency condition for the Chern-Simons 
coefficient k in the magnetic field. 



B. The improved photon propagator 

In Section [V] we shall investigate dynamical symmetry breaking by means of the 
Schwinger-Dyson gap equation for the fermion self-energy in the lowest Landau level. Ac- 
cording to the planer property of our model, the radiative effect of an improved photon 
propagator is saturated in the infrared region of its momentum, which will be shown in the 
following. 

Recall that the improved photon propagator D^ u , which is deduced from the effective 
action T[A] in Eq. ( |2.20| ), is given by 

D- u \p) = A- 1 u (p)-U IMU (p). (4.18) 
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The tensor form of D^ v becomes 



D^(p) = D e {p) 



PJIPu 

p\ 



+ D {p) ie^ + £ 



9^-^- + D ±(P) (I 
with a covariant gauge fixing £. The functions D e , Dj_, and D are determined as 

D e (p)g(p) = i 



(4.19) 



2 

n e (p) + ^n ± ( P ), 



fli(p)5(p)=4 n i(p). 

k - n (p) 



(4.20) 



D (p)Q(p) 



P 



where we have introduced the function which is related with a determinant of D 1 and 
is defined by 



Q{ P ) ■.= [i + n e (p)] { P 2 [i + n e (p)] +pin ± ( P )} 
+ [ K -n ( P )] 2 

= -e^-MetD" 1 ^). 



(4.21) 



The consistency condition, « = II o (0), means that the function Q{p) reduces to zero as 
p^ — > in the limit J± — > 0, where symmetries are restored in the tree- level Lagrangian. 
Therefore a photon becomes a massless propagating mode which is regarded as a Nambu- 
Goldstone (NG) mode due to the spontaneous breaking of Lorentz symmetry 0. 

Under a given magnetic field, we have shown that there exist three possible theories with 
the quantized CS coefficient \k\ = 0, a, and 2a, each of which has a massless photon. As far 
as we are concerned with the infrared behavior of the photon, it is enough to approximate 
the vacuum polarization functions with its values at zero momentum. From Eq. ( 4.12| ) and 
Eq. (|4.14 ) the parity-conserving parts U e and II_l are read as 



II e (0) = -6v / 2C(-l/2)«/, 
n ± (0) = 3v / 2C(-l/2) al, 



(4.22) 



which are irrespective of approaches to the symmetric limit J± — > 0. The improved photon 
propagator behaves, for all of the theories with \k\ =0, a, and 2a, like 



DM 



p 2 [l + cal] 



p2 



(p 



,2\2 



+ 



c [p±9iu-piPui 



with a positive constant c 



p 2 [l + cal] {2p 2 [l + cal]-cp 2 ± }' 
6 v / 2C(-l/2) ^ 1.76397. 



(4.23) 
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C. Expansions in infrared and ultraviolet regions of momentum 



In Section |VT] we shall investigate stability of vacua which are realized as possible solutions 
of SD gap equation, through the calculation of vacuum energy. For this purpose we need 
to know the asymptotic behavior of vacuum polarization functions both in infrared and 
ultraviolet regions of momentum. 

First we separate the functions into two parts according to p 2 < l~ 2 or p 2 > l~ 2 such as 

n,( P ) = n<( P ) e(i - (ipf) + n>(p) e((i P f - 1), (4.24) 

where i = e, _L, and o. The functions Hf (p) are estimated by expanding them with respect 
to (lp) 2 . In rigorous sense the expansion around p 2 = is possible only if the functions are 
analytic around p 2 = 0. Indeed at zero momentum the functions Uf (p) show singularities 
of delta functions when = £ n (n > 1) as read from Eq. (14.14]) . However, in our present 
study, since the parameter /i is nothing but the probe to pick up the spontaneous breaking 
of Lorentz symmetry and finally we take the limit J± — > 0, we only have to consider the case 
of \n\ < Si, which is free from the delta function singularities at the zero momentum. It is 
also the reason why we believe n o (0) is unaffected by higher-order corrections even when 
\fi\ = m. In the higher momentum p 2 > l~ 2 , the functions Hf(p) are expanded with respect 
to (lp)~ 2 by means of the asymptotic expansion. Since the magnetic field can be treated as 
a weak field in comparison with the momentum, the expansion with (lp)~ 2 matches the one 
with respect to the magnetic field. 

Now we shall provide a formal expansion for each case. In the lower momentum (strong 
field regime) we can formally expand exponential factors in Eq. ( |4.4| ) and Eq. ( |4.10| ) with 
respect to (lp) 2 . By using x := (lp) 2 , they are given by 

n<(p) = aJ [if+ii 1) x+L^ 2 + --" ) 

ITf (p) = al [Lf + L? J x + Lf x 2 + ■ ■ •] , (4.25) 
U<(p)=a[L^+L^x+ L^x 2 + --- 

All coefficients L&>, Lj , and (n = 0, 1, ■ • •) are dimensionless and show the delta 
function singularities in \fi\ = S n (n > 1), which means breakdown of analyticity around 
p 2 = 0. However, as far as we consider the symmetric limit J± — > 0, such singularities do 
not appear in any coefficient. Actually, our vacuum is realized in the symmetric limit as 
the lowest Landau level which is free from the delta function singularities. The first order 
coefficients are given by 

L^ = c, Lf = -°-, and = c, (4.26) 

where, for convenience sake, we have introduced c := k/ol which becomes c = 0, ±1, and ±2 
in the symmetric limit. 

In the higher momentum (weak field regime) the functions Hf (p) are systematically 
expanded as a power series of 1/x as shown in Appendix ||. They are formally given as an 
asymptotic series in the symmetric limit such as 
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#(°) + #w - + m-; _ 4 



r(l) 



(2) 



X 



(4.27) 



nT(p) 



a 



#(o) + I + ^(2) J_ + . . . 



According to the even or odd character under the parity, we find that coefficients H < f ,n+l \ 
H { l n+1 \ and vanish for all nonnegative integers n. We approximate the vacuum 

polarization functions as their asymptotic expansions up to 0(x~ l ) and the coefficients are 
given by (see Appendix [B|)[] 



H (o) 


= lim 




H (0) 


= lim s 




J±^0 




= 0, 



9(m - \n\)I(p; m) + 9(\p\ - m) 9 p 2 - 4(/i 2 - m 2 ) I(p; /i) 



m 



4 [p 2 - 2(/x 2 - m 2 )} 



PyP 2 ~ 4(/i 2 — m?) 



(4.28) 



where the function I(p; u) is given by Eq. ( |B15| ) in Appendix |B] and we use sgn(ft) = 
—sgn(fieB) , which obeys the Gauss law constraint [cf. Eq. (|3.20|) 1. 

We notice that Hf^ coincides with the coefficient calculated in the absence of the mag- 
netic field. This expansion in the weak field regime turns out to match the one around B = 
(or I — > oo). The coefficients are determined in the symmetric limit such as (see Appendices 
H and for details) 



7T 

2' 



H 



(0) 



0, and HP = 2c. 



(4.29) 



We will show in Section [V| that the coefficient , which is proportional to CS coefficient 
k, essentially contributes to a linear term of the magnetic field B in the shift of photon 
zero-point energy. 



V. SCHWINGER-DYSON GAP EQUATION 

In this section we study the Schwinger-Dyson gap equation for the fermion self-energy 
and show that dynamical breaking of both U(2N) and Lorentz symmetries are realized as one 



of nontrivial solutions to the SD equation. As confirmed in Section [[I IB], both condensates 
corresponding to the spontaneous breaking of U(2A) and Lorentz symmetries necessarily 
appear in such a complementary manner as to keep (ipip) — |("0 "0)1 a ^ a nonzero constant. 
It is a specific feature of (2 + l)-dimensional physics due to the kinematics of fermions in 
the magnetic field and this feature is also reflected in the structure of the SD equation. 

The fermion self-energy contains, as its dynamical variables, a scalar component m& and 
a 7° component /i^, which correspond to the dynamically induced fermion mass and chemical 



3 We use p as p = \fp* except for the usage of p as an argument of functions. 
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potential, respectively. These variables, and pd, are determined as solutions to the SD 
gap equation in the same self-consistent manner as the Hartree-Fock equation for energy 
gap in BCS theory. The analysis of condensates in Section [Tl rB| implies that it is convenient 
to use the proper combination J ± = p ± m instead of the explicit breaking parameters, p 
and m, in order to manifest the spontaneous breaking of symmetries. Subsequently, we are 
led to combinations; uj± := (pd ± m^l |18|]. This seems quite natural assignment once we 
recognize that — uj- corresponds to the energy gap of fermions on the lowest Landau 

level relative to the chemical potential (or Fermi energy) when pd > (pd < 0). In fact the 
coupled SD equations for rrid and fid can be reduced to two decoupled equations. One is the 
gap equation for uj + and the other is that for cj_ as shown later. This is another specific 
feature of (2 + l)-dimensional physics as that in condensates. 



A. Schwinger-Dyson equation in the ladder approximation 

We start our argument with construction of the SD equation in the ladder (bare vertex) 
approximation. In Section [TI| we have derived the 1/N leading effective theory in which a 
magnetic field couples to fermions as a background field. Our aim of gap equation analysis 
is to confirm realization of spontaneously magnetized vacuum supported by the condensate 
{ip^ip) 7^ 0. As in BCS theory, what sorts of excitations appear in a given system depends 
on the ground state which should be realized as one of solutions to the gap equation for 
self-energies of the excitations. 

In our present study the proper excitation is nothing but the fermion coupled with a 
magnetic field. Therefore we must construct the SD equation based on the bare fermion 
propagator S defined by Eq. ( |2.13| ) which should be recognized as the one for the proper 



excitation on the magnetized vacuum. As to the photon propagator, it has already been 
derived from the 1/N leading effective action (|2.20|) and has shown the massless photon 
behavior through the Gauss law as given by Eq. ( [L23| ). Thus the SD equation has the same 



contents as the one in QED3 with a magnetic field |17j apart from the fact that the fermion 
self-energy contains a 7 component other than a scalar component. 
The SD equation is given by the recurrent form of 

G{x, y) = S(x, y) - ie 2 f d 3 z dH S(x, z)^G{z, t)^G(t, y)D^{z - t), (5.1) 



(see Fig. ^ where G denotes a full fermion propagator which should be consistently deter- 
mined through the SD equation. We assume that the full propagator G also has the same 
form: 



G(x, y )=exp(^^(x- y)»A? (x + y)) G(x - y) (5.2) 



as S in Eq. (|2.14[) . We substitute the propagators into Eq. ( |5~T1) and perform the Fourier 



transform. After integrating out some of spatial coordinates, we obtain 
G(p) = S(v) - ( d 3 kd 2 q ± d 2 R ±e - iR ^ 

27T 5 J 



2 2 



) YG(k)YG^p°, P± ~ q f + 
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X 



D^(p° -k°,p ± -k ± + eBR 



(5.3) 



where := (0, R u R 2 ) and := (0, -R 2 , RJ. 

Following Ref. fll7| , we assume the strong magnetic field m {m d ) <C so that Eq. ( |5.3|) 
can be simplified owing to decoupling of the higher Landau levels S n (n > 1). The bare 
fermion propagator S(k) in Eq. ( |2.15|) can be decomposed into the Landau level poles 



5*(A;) = exp 
The function D n is given by 

D n (eB,k) = (m + k e7 °) 



E(-!) n 



D n (eB,k) 



n=0 



m? + 2\eB\n — k^ 



(5.4) 



leBl 



2A, 



3Ui 2- 



2A_ 



(5.5) 



where .£/"(:£) are the generalized Laguerre polynomials and A± denote the projection oper- 
ators to spin states: 



A, 



1 =p ij sgn(ei?) 



(5.6) 



We can easily see in Eq. ( |5.4Q that under the strong magnetic field, m <C the lowest 
Landau level dominates and all other higher Landau levels are negligible. Therefore, in SD 
equation ( |5.1| ), we approximate the Fourier transform of S with their lowest Landau level 
contributions, that is 



S{k) 



-i 2 k 2 



1 



m — 7° k, 



2A, . 



(5.7) 



We notice that fermions on the lowest Landau level essentially behave like (0+1) -dimensional 
objects. It is natural to write down the full fermion propagator also in the following (0 + 1)- 
dimensional form: 



G(k) 



9(k°) 2A 



+ ; 



(5.8) 



where g(k°) is a matrix commutative with A + . 

Substituting S and G into Eq. (|5.3|), we carry out q±,R± integration so that Eq. 



becomes simplified (0 + l)-dimensional form [17 

-,2 



rV) = "»-r(p°+A*) 



(2tt)< 



dk^gik^Dip - k°) 



(5.9) 



where the function D is defined by 

D(p°) :=-J <f P± e- l2 *>V 2 D 00 (p ,p ± ). 



(5.10) 
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Note that since we assume there is a magnetic field, the Gauss law (|1.1|) forces the effective 
CS coefficient at zero momentum to be zero. The photon in the magnetic field, therefore, 
behaves as a massless mode in infrared region no matter how we take the symmetric limit 
J± -> 0. According to Eq. fl423|), the one dimensional photon propagator D is represented 
by the integral exponential function which has a logarithmic behavior in the infrared region 
of momentum such as 

^LD(p) n -?±)n\lp\, (5.11) 

where «o := al/N(l + cod) and p := — ip°. The constant c denotes parity even vacuum 
polarization effect, n e (0) = cal, and is given by c = -6v^C(-l/2) - 1-76397. Thus only 
the infrared momentum region is relevant in Eq. (|5.10|) . 



B. Rearrangement of dynamical variables 



The matrix function g is written in an SO (2) invariant form; 

g-\p°) = B(p°) + A(p°) p° - 7° \B(p°) + A{p°) p° + te sgn(p°; 



(5.12) 



where the functions A, B, A, and B are even functions of p°. The self-energy part is divided 
into scalar and 7 components as the functions B and B, respectively. Since p° is SO(2) 
invariant in itself, A(p°) p° should be involved in g~ x as an odd function part of the scalar 
component. We can set the functions A(p°) and B{jP) to be positive definite without loss 
of generality. Following the physical implication mentioned in the beginning, we rearrange 
the functions A, B, A, and B as 



A±(p°) := A(p°) ± A(p°) 
B±(p°) :=B(p°)±B(p ) 



(5.13) 



If we recognize that the functions ±B±/A± are identified with the relative energy of the 
proper excitation to the Fermi energy, it seems natural to assume that B±/A± have a 
definite sign irrespective of its argument. Then we can perform the Wick rotation k° = ik in 
the SD equation without any ambiguity^ Actually g(k°) is decomposed into two propagators 
which correspond to the proper excitations in the lowest Landau level; 



1-7° 
2 

1 + 7° 



1 



B + (k°) + A + (k°) k° + iesgn(k°) 
1 

B_(k°) + A_(k°) k° + iesgn(k c 



(5.14) 



whose poles are located on a complex k° plane at 



We use k, p instead of k, p only in this section to simplify the mathematical formulas. 
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\B±/A±\ + ie for B±/A± > 
\B±/A±\ -ie for B±/A± < 



(5.15) 



Thus we can perform the Wick rotation without any residual contributions from poles. 

Substituting Eq. Q5.12|) into the SD equation ( |5.9| ) and performing the Wick rotation 
k° = ik, p° = ip, we obtain two sets of coupled integral equations |18 



B±(p) - J± 
p[l-A ± (p)] 



(2tt) 3 

„2 



(2 



7T 



dk 



dk 



B ± {k) 



A 2 ± (k) k 2 + Bl{k) 

A ± {k)k 
A 2 ± {k) k 2 + B 2 ± {k) 



D(p-k), 
D(p-k), 



(5.16) 
(5.17) 



which are divided into the coupled equations for (A + , B + ) and those for (A-, BJ) by virtue 
of the specific feature of (2 + l)-dimensions [cf. Eq. ( |5.14|) 1. It is obvious that the above 
equations have trivial solutions B±(p) = in the symmetric limit J± — > 0. 

Firstly we determine A±(p) for the trivial solutions, B±(p) = 0, as well as for the 
nontrivial solutions. We differentiate both sides of Eq. ( |5.17| ) with respect to p and obtain 



l-A ± (p)-pA' ± (p) = -^£ 



dk 



A ± (k)k 



A 2 ± {k) k 2 + Bl{k) p-k' 



(5.18) 



where we have used the asymptotic form of D in Eq. (|5.11 ). According to massless photon 
behavior in the infrared region, the integral in Eq. ( |5.18|) shows infrared divergence if we 
naively set B±(p) = 0. Besides a spurious photon mass and a gauge fixing parameter do not 
play any role to regularize this infrared divergence because of the logarithmic behavior of 
D(p). We therefore need to leave B± in the integrand as a cut-off to regularize the infrared 
divergence even for the trivial solutions. 

Let us suppose that A±(p) are constant in almost all of the momentum regions except 
for p«0, while pA'±(p) are negligible even in p w 0. For the large momentum p the integral 
in Eq. (|5.18|) vanishes and we get A±(oo) = 1. For p w we approximate the fermion 
propagator as 

(5.19) 



A 2 ± (k) k 2 + Bl(k) k 2 + Bi(oy 

which obeys our assumption. Thus Eq. (|5.18|) reduces to 



At(O) 



Til 



k 2 + Bl(0) 



(5.20) 



and we obtain [[H]] 



A ± (0) 



1 + 



Of, 



\IB A 



(5.21) 



Although A±(0) reduce to zero for the trivial solutions B±(p) = 0, A±(p) ~ 1 are satisfied 
except for the momentum around p ~ and are consistent with the requirement of Ward- 
Takahashi identity under the ladder (bare vertex) approximation. 
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Now we solve Eq. (|5.16|) and find the nontrivial solutions B±(0) 7^ in the symmetric 
limit J± — > 0. Since the integral in Eq. ( |5.16|) is dominated at the infrared region by the 
logarithmic behavior of the photon propagator in Eq. ( |5.11| ), the approximation in Eq. ( |5.19| ) 
is also valid and Eq. (|5.16|) is simplified as 



B±(0) 



dk 



B ± {k) 



k 2 + Bl(0) 



\n\-lk\. 



(5.22) 



If the above gap equation has nontrivial solutions, they should satisfy |Zi?±(0)| C 1 so as to 
be consistent with the lowest Landau level dominance. The smallness of the solutions leads 
to the result that the dominant contribution comes from the infrared region, k ~ 0, in the 
integral in Eq. ( |5.22j ). Thus we can replace B±{k) with B±(0) and obtain the gap equation: 



a 

u± = 

7T 



00 o ~\- UJ± 



ln I — s\ 



(5.23) 



where we introduce dimensionless variables u± 



solutions ±u s given by u s = —ao lnu s , as well as the trivial one [17 



IB±(0). This equation has the nontrivial 
Note that the nontrivial 



solution lu s satisfies the condition u s < 1, since ao < 1 for any e and N. If u± appear as 
one of the nontrivial solutions ±u s , A±(0) are determined as 



A ± (0) 



1 - 



Into* 



-1 



1. 



(5.24) 



which is consistent with the Ward-Takahashi identity, A±(p) = 1, due to the smallness of 
the nontrivial solution a;,. 



C. Classification of the nontrivial solutions 

The obtained gap equations are the same for all three theories assigned to \k\ = 0, a, 
2a. On the other hand, in order to maintain the self-consistent magnetic field for each k, 
the dynamical mass and the chemical potential \id have to obey the same relation as 
the explicit breaking parameters m and \i satisfy in Eq. ( |3.2U| ) (see Fig. |2|). In fact u± do 
not necessarily choose the nontrivial solutions ±u s . Consistent solutions are automatically 
assigned to each k such as 

!(UJ S , -U s ) , K = 

(u s , 0) or (0, -u s ) ,\k\ = a , (5.25) 
(u s , u s ) , |«| = 2a 

where the replacement of u s with —uj s provides another set of nontrivial solutions. It is 
intriguing that the solutions for \k\ — a appear as combinations of a nontrivial solution and 
a trivial one. The dynamical variables, and fid, which are defined by := (u + — u-)/2l 
and fid := (u + + u-)/2l, respectively, are determined as 

f (u s /l, 0) for k = 
(m dl \fi d \) = I (u s /2l, u s /2l) for \k\ = a . (5.26) 
[ (0, u s /l) for \k\ = 2a 
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For k = 0, fermions acquire their dynamical mass, while the vacuum becomes empty 
(fid = 0). The fermion self-energy changes entirely into the dynamical mass and reproduces 
the result in Ref. |T7|. For \k\ = 2a, the self-energy is used to occupy the lowest Landau 



level fully with fermions. Then there is no mass generation and it supports dynamically 
the result in Ref. |J. Noteworthy case is \k\ = a where the self-energy is shared by 
and fid- The vacuum is realized as the lowest Landau level half-filled by massive fermions. 
This solution is a new one which spontaneously breaks the U(2iV) symmetry as well as the 
Lorentz symmetry at the same time ||I8|. 



VI. SPONTANEOUS MAGNETIZATION AND FERMION MASS GENERATION 

It still remains unknown whether or not the nontrivial solutions to SD gap equation 
as well as the self-consistent magnetic field are energetically more favorable than a trivial 
solution without any magnetic field. In order to solve the above problem, we have to 
investigate the effective potential which is obtained as the vacuum energy shift owing to the 
spontaneous breaking of Lorentz symmetry and/or XJ(2N) symmetry. Since the solutions, 
irid and fid, are connected with the self-consistent magnetic field B, the effective potential 
depends only on B and is given as a function V(B). It is composed of four parts including 
the Maxwell energy; 

V(B) = V CJT (B) + V F (B) + V P (B) + ^ . (6.1) 



Vcjt denotes the CJT potential |20] which gives the energy difference between a nontrivial 



vacuum and a trivial one under the presence of the magnetic field. Vf (Vp) corresponds to 
the shift of fermion (photon) zero-point energy due to the magnetic field in the symmetric 
limit. 

In the following we derive the potentials, Vcjt, Vf, and Vp, separately and investigate 
stability of vacuum in the large N limit for each value of k. We employ strong coupling 



expansion with respect to 1/al ~ ^B/e 3 in order to estimate the lowest order contribution 
to Vp based on the expansion of vacuum polarization functions in two regimes of momentum, 



which is shown in Section IV 



A. Vacuum energy shift due to fermions: Vf 

The vacuum energy in the large N limit is given by the first term of effective action in 
Eq. (|2.20| ). It depends on the explicit breaking parameters (m, fi) as well as the magnetic 



field B such as 

VLEi(m,fi;B) := iNTrLnS' 1 , (6.2) 

where Q denotes the three-dimensional space-time volume and S -1 is the inverse of a bare 
fermion propagator. It is nothing but the zero-point energy for fermions coupled with the 
magnetic field B. The shift of zero-point energy for fermions due to the magnetic field is 
determined in the symmetric limit as 
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V F (B) := lim [E x (m, //; B) - £7i(m, //; S -»• 0)] . 

J±— > o 



(6.3) 



It is convenient to use the condensate Jo to derive the finite density part of E% (/i 7^ 0) 
|33|. From the definition ( |6.2|) , we notice 



dEi(m, fi; B) 



-n-HNTi[-f°s} 

-Jo(m,n), 



(6.4) 



where Jo is the condensate calculated in Section [TT1] and given by Eq. ( |3.13| ). Integrating 
both sides over /x, we obtain another form of E±; 



Ei(m, //; B) = Ei(m, 0; B) — / dn'Ja(m, //; B). 



(6.5) 



As to the computation of zero density part, the proper time method is used. We see from 
the definition (16.21) that 



VL E x (m, 0; B) = iN Tr Ln [m - iV] 
= — Tr Ln [m 2 + V 2 ] 



(6.6) 



where V := 7 /i D A1 [A ext ] and we have used the relation 7 5 V7 5 = —V. The symbol Ln denotes 
a logarithmic function whose argument is given by operators, and Tr means a trace over 
space-time coordinates as well as over spinor indices. The TrLn term is rewritten with an 
integration over the proper time; 



Tr Ln [m 2 + V 2 

Eqs. (|2.13| - |2~T5|) lead to the following relation 

d 3 k 



(l'\r I — tr (x\ e 

1/A 2 s 



-is [m 2 +V 2 



\x). 



(6.7) 



(x e 



-is [m 2 +V 2 



\x) 



(2vr) : 



exp 



is m 2 - k 2 Q + tan(eEg) k 2 
V eBs 



1 + 7 V tan(e£s) . (6.8) 



Applying this relation to Eqs. ( |6.6[ - |6.7| ), we obtain, after Wick rotation and Gaussian inte- 
gration, 



£?i(m,0;5) 



N 



4vr 3/2 

N 



\eB\ 



ds 



1/A 2 s 3 / 2 
N 



e~ sm coth(|e£|s) 



— --Am 2 

2tt 3 / 2 2tt 



efil 3 / 2 



m/ + 2\/2C 



1 (m/) s 



+ 1 



(6.9) 



where the function ((z, q) denotes the generalized Riemann zeta function which has appeared 
in Section |T|. The finite density part is calculated from the condensate J in Eq. ( p,13j ) as 



f dp,'Jo{m,v';B) = ^-\eBf 2 
Jo Ztt 



- m)9(\n\ - m) + 2 £ l(\fi\ - S n )6(\fi\ - S n 

n=l 



(6.10) 
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Thus the explicit formula for E\ is given by 



Eiijn, //; B) 



N k n N . ,o /9 

— — Am 2 eB 3/2 

2vr 3 / 2 2vr 1 1 



m 



/ + 2V/2C 



1 {mlf 



+ - m)0(M - m) + 2 ^ /(H - - £*) 



ra=l 



(6.11) 



The fermion vacuum energy in the absence of the magnetic field is obtained by taking the 
weak field limit in Eq. (|6.11|). We find 



E x {m,n;B -> 0) 



N a 2 W 3 

— —Am H m , 

2tt 3 / 2 3tt ' 



(6.12) 



where we have used the asymptotic formula for the zeta function p5| , 

1 



l)q 



2-1 



2q 



oo 



(6.13) 



We notice that the ultraviolet divergence which appears in B ^ can be completely canceled 
out by another ultraviolet divergence in B = 0. 

After taking the symmetric limit, Vp is determined as an ultraviolet finite function of B, 
that is, 



V F (B) = -^-\eB\ 3 / 2 4V2((-l/2), 



(6.14) 



which is irrespective of the values of k, or approaches to the symmetric limit. 



B. Vacuum energy shift due to the nontrivial solution: Vcjt 

The Schwinger-Dyson gap equation is derived from the effective action for composite 
operators, ipip and It is given by a functional V[G] of the full fermion propagator G 

which is defined by 



G{x,y)5 ij :=i{0\Tijj i {x)i; j {y)\0). 



(6.15) 



The equation which leads to stationary points of T[G] is equivalent to the SD gap equation. 
The effective action T[G], or Cornwall- Jackiw-Tomboulis action, is given by pO 



NT[G] 



-iN Tr 



LnG- 1 + S^G - ll +NT 2 [G], 



(6.16) 



where T 2 [G] denotes all of the 2PI (two-particle irreducible) bubble diagrams composed of 
the full propagator G. The large contributions to T 2 [G] are given by a gauge invariant 
two-loop diagram (see Fig. |]); 



ivr 2 [G] 



Ne 2 



Tr [GVGY] D 



(6.17) 
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where the photon propagator D^ v involves the large N screening of fermions and has already 
been given in Section [V. We see that a stationary point for T[G] actually provides the SD 
gap equation 



5T[G] 
5iG 



G 



-1 



S 



-1 



sr 2 [G] 

5iG 



0. 



(6.18) 



which has been solved in Section [V| and shows the nontrivial solutions given by Eq. (]5.26|) 
as well trivial one. 

What we want to know is stability of the vacuum which is realized as a nontrivial solution 
(m d , Hd) for each k. To see this we only have to investigate which is energetically more 
favorable between the nontrivial solution and the trivial one. The vacuum energy due to the 
dynamical solutions (m d , fi d ) is provided by the CJT action as its value on the stationary 
point; 



fi%,wB) :=-iVr[G so i] 



ziVTr 



1 



L n G S oi + - (S 1 — G en { ) G 



sol 



'sol 



(6.19) 



where G so i means the solutions to the SD equation ( 6.18|) and we have eliminated the 2PI 
contribution by means of Eq. ( |6.18| ). Thus the vacuum energy shift due to the nontrivial 
solution is determined as 



Vcjt(S) := hm [E{rn d , m B) - E(0, 0; B)} . 



(6.20) 



To estimate Vcjt from the solutions obtained in Section [V], we adopt the following procedure 
as an approximation, that is the replacement 



lim G sol ~ 



S\ 



(6.21) 



where S denotes the bare fermion propagator defined by Eq. ( |2.13|) . All of the Landau levels 
contribute to Vcjt- However, the lowest Landau level dominates due to smallness of the 
solutions, m d l 1 and [ij "C 1. Under the replacement ( |6.21| ) we see 



lim iN TrLnG^ = nE 1 (m d , f i d ; B) 



(6.22) 



where the energy function E\ is the same as the shift of fermion zero-point energy. We also 
find that the second term in Eq. ( |6.19|) becomes 



hm — Tr Ln 

j±^o 2 



(/i d 7° - m d ) G so \ = ^ J (m d , fi d ) - m d J(m d , fi d )] , (6.23) 



where the functions Jo and J are the same condensates calculated in Section |TJ except for 
the replacement (m,n) — > (m d ,n d ). Using Eqs. fl3.12| , |3.13"D and Eq. (|6.11| ) we obtain 



lim E(m d} n d ;B) 
j±— > o 



-— \eB 3/2 

47T 1 1 



m d l + l(\/j, d \ - m d ) 9{\ji d \ - m d ) 



-V2(m d iy c (l^ 



AV2C 



1 {m d lf 
2' 2 



(6.24) 
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where only the lowest Landau level contributes to the finite density part due to smallness 
of the solutions, while the zero density parts appear as a summation over all the Landau 
levels. Notice that the ultraviolet divergences in E\ and J are completely canceled out in 
the whole of E. 

Thus Vcjt is determined in such a form which depends on the nontrivial solutions in Eq. 
( pg ) as 



VcMB) = -^\eB\^ 
Air 



lmax{m d , \fjL d \} + 0((m d iy 



(6.25) 



which shows that the vacuum energy is shifted to be negative. Therefore it is confirmed that 
dynamically generated solutions are energetically favorable for each k in the self-consistent 
magnetic field B. 

C. Vacuum energy shift due to photon: Vp 



Vp is derived from the effective action (|2.20| ) as the next-to-leading-order contribution 
in 1/N, that is, 



where we have introduced D^ip) which is defined as the photon propagator D^ u {p) in 
B — > 0, and used matrix notation for space-time indices. The symbol tr thereby means the 
trace for space-time indices. We notice that Vp is the shift of photon zero-point energy due 
to the magnetic field.0 

In order to estimate the integral in Eq. (|6.26|) , we have to expand Vp into a power series 



with respect to the magnetic field B. However we have supposed the strong magnetic field 
which supports the lowest Landau level dominance and therefore guarantees the analysis of 
the SD gap equation in Section [V]. In fact we are led to the small solutions, m d (\fi d \) <C 
in comparison with the magnetic field. Hence we need another scale which can be taken to 
be larger than the magnetic field and provides the small parameter for the expansion. 

We have indeed such a scale, that is, the gauge coupling a to which the Chern-Simons 
coefficient k is related through the Gauss law. The parameter r := 1/al can be small by 
choosing a sufficiently large a with the magnetic length I fixed finite. This results in a 
hierarchy of three scales, m d (\fi d \) <C / _1 <C a, which should be achieved for our present 
study and the solutions {m d , fi d ) actually allow this situation. Thus we recognize that the 
expansion of Vp(B) is a kind of strong coupling expansion with respect to 1/a. It is enough 
to pick up and estimate the lowest order contribution of the expansion with respect to r in 
order to investigate stability of vacuum. 



8 Hereafter we use the propagator D (D) where the symmetric limit J± ^ has been taken. 
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To do the above procedure it is necessary at first to expand an integrand of Eq. (|6.26 ) in 
such a consistent way that matches the expansion with respect to 1/od. In the following we 
apply an expansion of D around D and try to estimate the lowest order contribution to Vp. 
The analysis is based on the expansion of vacuum polarization in two regimes of momentum 
in Section [TV[ 

The propagators D and D are given by 



D-l{ P ) := A;J(p)-n^(p), 

D-J(p) :=A;»-fV(p), 



(6.27) 



where (p) denotes the vacuum polarization in the limit B — > and is expressed as0 



IV (p) = (PuPv ~ p 2 9^) n(p), 



(6.28) 



which is written in the gauge invariant form even in the presence of the finite density fi ^ 0. 
The function H(p) is determined as 



n(p) 



lim — 

n\ a 



9(m - \(i\)I(p; m) + -m)9 p 2 - 4(/i 2 - m 2 ) I(p; //) 



2/ p 



(6.29) 



with the function J(p; m) given by Eq. (|B15| ). From Eq. ( |6.27 ) D is rewritten in the form of 
expansion around D such as 



6U^(p) :=U^(p)-n^(p). 
Thus we obtain the following expansion 



(6.30) 



tr Ln 



D(p)D-\p)] =trLn[l-D(p)<5n(p) 



OO 1 

£-tr \D(p) 5U(p) 



n=l 



n 



(6.31) 



Each term in the above summation is estimated through the decomposition of a rank 2 3x3 
matrix D 511 into a 2 x 2 unit matrix and Pauli matrices. It is expressed as 



tr 



D(p)5U(p) 



9 [ ^ ] r [R( P )r 2m [s\p)r 

m=0 



[Q(p)f 



(6.32) 



where [n/2] means the integer part of n/2 and the functions Q, R, and S = (Si, S2, S3) are 
determined as 



9 The parity-violating part vanishes in the limit B — > as shown in Eq. ( (4.10| ) due to our setting 
of fermions. 
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q(p) = p 2 i + n(p) 



i?(p) = -p 2 i + n(p) 



2 

5n e (p) + ^5n ± ( P ) 



2p 2 



+ k5II (p), 



Si(p) 



i + n(p)] ^tfn ± (p), 



(6.33) 



pi 



5 2 (p) = g K OT ± (j9), 

s 3 (p) = -i P { [i + n(p)l <$n (p) + k 



2 

OT e ( P ) + |^<m ± (p) 



The deviations 8H e , SU±, and «5I1 are the counter parts for IT e , IT^, and n o in Eq. fl6.30| ), 
respectively. 

In Section [TV] we have expanded the vacuum polarization functions in each regime of 
momentum. The energy shift Vp is divided into two parts according to two regimes of the 
loop momentum. After the Wick rotation, it is given by 



oo [n/2] 

v P (B) = E E 



n ° 2m [V< m (B) + V> m (B)] 



V< m (B) :-- 



n=l m=0 
-1 



n 



47T 2 

-1 

47T 2 Jo 



d6 sin 6 I dpp 2 
Jo 

7T /*00 

d#sin0 / dpp 

Ji- 1 



i- 1 - [^(p)]"" 2 " 1 [S 2 (p)] m 



[Q(p)] n 

[i?(p)]- 2m [S 2 ( P )] m 



(6.34) 



where 9 denotes the polar angle for the Euclidean momentum We apply the (asymp- 
totic) expansion for the vacuum polarization functions in the strong (weak) field regime of 
momentum to V< m (V> m ). 

In the strong field regime (p < / _1 ), we notice from the results of Eq. (|4.26|) that 



<x 



71 



- + y/xL e {x) 



5H± = al L±(x), 5H = a L (x), 



(6.35) 



where L e , L±, and L Q are expressed as analytic functions. Substituting the above results 
into Eq. ( |6.33| ), Q, R, and S 2 are written in such a factorized form as 



Q = a 2 



7T 

- + ry/x ) + c" 



R = a 2 f\\fx\ r 



e 2 4 
o = a 



g(y/x\ r) 



(6.36) 



where analytic functions f(u; r) and g(u; r) depend on the small parameter r := 1/al. Then 
V< becomes 



2vr 2 / 3 



1 ^ d ,/ 1 du « 2 [/(«;r)r 2 ™W«;r)r 



Jo 



[ru + f ) 2 + c 2 



(6.37) 



We notice that a factor — 1 /27r 2 Z 3 supplies r 3 . If an integral over u diverges by going to 
r = 0, it lowers the power of r to be less than three. Since in the limit r — > both of / and 
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g are analytic around u = and denominators of integrands become constants, the integrals 
in Eq. (|6 . 3 7| ) converge and the lowest order is determined as 0(r 3 ) which is the same order 
as Vp and Vcjt- 

In the weak field regime (p > / _1 ) we find from the results of Eq. (|4.29| ) that 



<5IL 



al hJ 1 - 



x 



5 / 2 e \x, ' 



5U 



01 Hi a 



X 



5/2 



X 



>)U„ =-hJ- 



x 



x 



(6.38) 



where H e , H±, and H Q are also expressed as analytic functions. Substituting the above 
results into Eqs. (|6.33|), Q, R, and S 2 are written in such a factorized form as 



Q = a 2 x 



7T 1 

2Vi 



| 2 + ^ 

X 



R 



or 



x 



X 



- G 7= : 

X \ v X 



(6.39) 



V> 



by using analytic functions F(u;r) and G(u;r). Then we obtain the following formula for 

u 4(n-l) r )]»-2m [ G ( u;7 .)]« 



U 



2 in 



> + f M) 2 + C 2 M 2 



(6.40) 



The convergence at r = is not so trivial this time because of the factor u i( - n ~^ / u 2m . 
However G(u; r) supplies another factor u 2 when r = 0, which can be seen in the explicit 
form of it. Therefore we obtain an identity G(u; 0) = u 2 G(u) with the analytic function 
G(u). Inserting the above identity into Eq. ( |6.40|) , we find in the limit of r — ► 



V> m (B) oc -± £d V f'duu^-V [F(u;0)} n ~ 2m [G(u) 



(6.41) 



The lowest order contribution would appear in the above integrals if they show divergence. 
The integral shows divergence only when n = 1, namely, V^. In order to extract the lowest 
order contribution, it is therefore sufficient to estimate the coefficient of the most singular 



>. 



part of the integral in V w . 



1 



F(u] r) 



,i ,i 

_ *h dr] Jo dU \r + luf + c 2 u r 



(6.42) 



From the definition in Eq. (|6.39| ) the function F(u; r) is given by 



F(u; r) = cH a (u 



H e (u 2 ) + ^H ± (u 2 ) 



(6.43) 



where we see that the second term supplies a factor u 2 at r = so that it never contributes 
to the lowest order. Thus it is clarified that the lowest order contribution is determined only 
by H o (0) such as 



2tt 2 / 3 
a 



1 /2c 

-H o (0) arctan I ) + const. + O(lnr) 

r \2r + 7r,' 



2tt 2 / 2 



H o (0) arctan 



2c 



7T 



(6.44) 
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This means that even if we take all of the large N contributions in vacuum polarization into 
account, the lowest order contribution to Vp is determined only by H o (0). Therefore the 
above result is exact in the large N limit. The lowest coefficient H o (0) has already estimated 
as H o (0) = 2c which leads the coefficient of \eB\ to be negative definite. 



Thus the vacuum energy shift Vp up to 0(B 3/2 ) is obtained as in Eq. ( ggg ). Specif- 
ically, it is given by 

V P (B) = -J4 \eB\ arctan [ ^ ) + 0{B 3 / 2 ). (6.45) 
7r 2 y na J 

It is obvious that the negative linear term for \eB\ is owing to the absence of topological 
photon mass in B ^ as well as the presence of it in B = as pointed out by Hosotani || . 
If we set c = ±2, the above result completely matches Ref. ||. 



D. Stability of vacuum 

Now we are in a position to discuss the stability of vacuum, or a possibility of spontaneous 



magnetization. Recall that the entire effect of vacuum energy shift is given by Eq. (|6.1| ). We 
notice that the negative linear term for B appears only in theories with nonzero K, in which 
the magnetic field is possibly supported by the charge condensation e(ip^ip). Therefore, in 
the theories with \k\ = a and 2a, the potential has its stable stationary point at B ^ 0, so 
the spontaneous magnetization occurs and Lorentz symmetry is dynamically broken. 

In the theory with \k\ = 2a, the SD gap equation tells us that the fermion self-energy is 
exhausted to induce a charge density on vacuum so that there is no generation of fermion 
mass. Then the vacuum appears as a fully-filled magnetized vacuum [El. Alternatively, in 
the theory with \k\ = a, the SD gap equation shows not only the charge condensation but 
also the dynamical fermion mass, whose stability under the background magnetic field is 
supported by the CJT potential. Moreover the entire energy shift shows that the background 
field itself is also induced dynamically. Thus we are led to the half-filled vacuum in which 
both of the U(2iV) symmetry and the Lorentz symmetry are spontaneously broken |[L8|| . 

As to the theory with k = 0, that is QED 3 with N four-component fermions, the situation 
becomes rather subtle because there is no Chern-Simons term, or the negative linear term in 
B. The lowest order term in Vp is provided 

£3/2 

term, so all of higher-order corrections 
in the expansion of contribute to its coefficient. However, in the entire energy shift 
V, the B 3 ! 2 term is saturated by that of Vp with a factor N. Therefore at least in the 
large N limit the entire potential has a positive B 3 ! 2 term and results in no spontaneous 
magnetization. 

We can see this more explicitly in the formula for Vp. Through an identity Tr Ln = In Det, 
it is rewritten as 



Vp(B) = lim - [-^- ln 



[l + ne^ni + n^ + ^i/^iij 



[i + n(p)p 



(6.46) 



in QED 3 [cf. Eq. (|4.21 )]. Since in QED 3 a photon field is always massless whether a magnetic 



field is turned on or off, there is no essential difference due to the magnetic field in the weak 



32 



field regime of momentum.^ We therefore cut off the integral at p = I 1 and approximate 
the vacuum polarization functions as their values at p = 0. Then we obtain 



V P (B) 



4tt 2 / 3 



dr] / du u 2 In 



[r + c][r + c — cr) 2 /2] 



[r + ir/2u} 2 



(6.47) 



We see that the integrand becomes negative around u = due to the \ jp behavior of H(p) 
so that Vp is led to be negative. The integrals can be performed analytically and Vp is 
obtained as 



V P (B) 



1 



\eBf*K(— 
12tt 2 ' 1 \wa 



eB\ 



(6.48) 



with the function K(x) defined by 
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+ In 



'x + b/2 s 
. x + 1 , 



^larccothj^). 



(6.49) 



where we have introduced b := 2c/tt ~ 1.12298. 

K(x) shows monotonically increasing behavior and becomes negative definite in x > for 
b < 4/3, that is K(0) < K(x) < as shown in Fig. |. Hence, if an inequality \K(0)\ < 2tccN 
is satisfied, then the B 3 ^ 2 term in the entire energy shift V(B) becomes positive semi-definite 
and shows monotonically increasing behavior which leads us to a non- magnetized vacuum. 
If we suppose large N, the above inequality is trivially satisfied, which only means that the 

£3/2 

term is dominated by the large N contribution of Vp. However the parameter c (or b), 
which is determined by the kinematics of fermion in the magnetic field, satisfies 



mo)\ 

2ttc 



2vrc 



— + 2V2 arccoth\/2 + \n(2c 2 /n 2 ) 



0.0572877. 



(6.50) 



Therefore, even if iV = 1, the spontaneous magnetization or equivalently the spontaneous 
Lorentz symmetry breaking does not occur in QED3 with N four- component fermions. 



VII. CONCLUSION 



In this paper we have investigated dynamical symmetry breaking in Chern-Simons QED 3 
associated with the realization of finite density vacua (or occupied lowest Landau levels). 

Through the explicit estimation of condensates (ipip) and (ip'ip), we have clarified that 
in the presence of the Chern-Simons interaction the magnetic field is necessarily connected 



This situation changes drastically in Chern-Simons QED3. If B = 0, a photon has a large 
topological mass of the order of a (> Z -1 ) without any screening. Therefore the photon behaves 
as a massive field even in the weak field regime of momentum (p > C 1 ) when B = 0. 
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to the charge condensation e{ip^ip) through the Gauss law in such a way as to restrict the 
possible values of the CS coefficient k to ±Ne 2 /2ir, ±Ne 2 /4ir } and fl8|,p8U29f . In other 
words, the magnetic field can penetrate the system only if the photon effectively becomes 
massless |§. We also have found that both condensates (ipip) and {ip^ip) always appear so 
as to complement each other and to keep the combination (tpip) — constant. The 

symmetry breaking patterns, or the vacuum configurations, are displayed in Table. I. 

In the theory with \k\ = Ne 2 /2ir, the vacuum is realized as a fully-filled lowest Landau 
level \(v + ,vS) = (1,0), or (0,1)] which is a singlet of the flavor U(2iV) group. Only the 
Lorentz symmetry is spontaneously broken by the induced magnetic field together with 
charge condensation, as pointed out in Refs. P|J. Alternatively, in the theory with k — 0, 
or QED3 with N four-component fermions, the vacuum is realized as an empty lowest 
Landau level [u± = 0] even in the presence of the magnetic field. The U(2iV) symmetry is 
spontaneously broken, while the magnetic field is not accompanied by charge condensation 
[TD| , I7| . As to the theory with |k| = iVe 2 /47r, both condensates coexist in the half-filled 



lowest Landau level = (1/2,0), or (0,1/2)]. This vacuum allows the spontaneous 



breaking of both symmetries, that is, the U(2iV) symmetry and the Lorentz symmetry |I 

The above situation suggests that in the fermion self-energy the scalar component rrid 
and the 7 component fid also complement each other. We have also attempted to find 
the dynamical solutions which correspond to the possible configurations of the above three 
classes through the analysis of the Schwinger- Dyson gap equation. 

In the theory with \k\ = Ne 2 /2n, the realization of a fully-filled vacuum is supported 
by the solution (m^, \fid\) — (0,^/0- The fermion self-energy has only a 7 component, 
the originally massless fermions remain massless and there is no generation of fermion mass. 
This solution reproduces and verifies the result of Ref. ||, that is, vanishing fermion mass 
and broken Lorentz symmetry. Thus, the magnetic field does not necessarily lead to mass 
generation if the vacuum is fully-filled by fermions. This is in sharp contrast to the results of 
Ref. [17 1, which are based on the empty vacuum and correspond to the theory with k = 0. In 
this case the solution to the gap equation has (rrid, \fid\) = (u s /l, 0). The fermion self-energy 
is saturated by the dynamically generated fermion mass, which then results in an empty 
vacuum. 

In the case \k\ = Ne 2 /4n, we have found a novel solution (m^, \fMd\) = (u s /2l,u s /2l) 
which causes spontaneous magnetization as well as the dynamical generation of fermion 
mass | 18fl . The scalar component rrid and the 7 component fid are comparable to each other. 
The fermion self-energy is shared, half filling the lowest Landau level and half supplying 
fermions their dynamical mass. In the resultant half-filled vacuum the U(2iV) symmetry 
and the Lorentz symmetry are simultaneously broken through the dynamically generated 
fermion mass and the induced magnetic field [0 . 

We have investigated the question of vacuum stability for each class through the explicit 
calculation of the energy shift V(B) due to the magnetic field B. We find that a negative 
linear term in B always appears for \k\ — Ne 2 /27r and iVe 2 /47r, leading to a nonzero magnetic 
field at the stationary point. Moreover, the estimated coefficient for the above term is exact 
in the large N limit and matches the former result in Ref. ||. 

For k = 0, or QED3 with N four-component fermions, the lowest term in the energy 
shift behaves as B 3 ^ 2 . The fermions and the photon make opposite contributions. Since the 
energy shift due to the fermions is the leading one at large N and is positive, the entire 
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energy shift becomes positive at least as iV — > oo. However, our rough estimate has shown 
that even if AT = 1 the entire shift is positive due to the kinematics of fermions in the 
magnetic field. Therefore we can conclude that spontaneous magnetization does not occur 
in QED 3 with N four-component fermions.^] 

It is an open problem whether or not the magnetized vacua found at zero temperature are 
maintained also at finite temperature, although the Lorentz symmetry is explicitly broken by 
the heat bath. The stability of the magnetized vacua at finite temperature was confirmed 
in Ref. [37] in the same context as Ref. ||. On the other hand, Ref. |38| found that in 



2 + 1 dimensions and in the external magnetic field the condensate {ipip) is unstable at finite 
temperature T due to nonanalyticity at the origin in the (m, T) plane. Thus one needs to take 
into account the possibility of dynamical generation of fermion mass as well as spontaneous 
magnetization also in the finite temperature. Another problem is to extend our study to 
non-Abelian gauge theories, namely, Chern-Simons QCD 3 , whose intimate relationship with 
the frustrated Heisenberg antiferromagnets is studied in Ref. |3Df . 

The gauge sector of Chern-Simons QED 3 is dual to the Abelian Chern-Simons-Higgs 
theory, where the kinetic term for the Higgs field provides the gauge boson mass term 
through unitary gauge fixing |4(|. On the other hand, the (2 + l)-dimensional Thirring 
model can be reformulated as an Abelian Higgs theory with the realization of a massive 
composite gauge boson through the introduction of a spurious Higgs boson f¥I|j . It may, 
therefore, also be interesting to study another infrared sensitive theory; 

C! = ^ [id„ + eAJi ip - 7^A^ - ^e fll/p A^A p , (7.1) 

whose dynamical gauge boson might play an important role in planar condensed matter 
systems. 
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Some authors [34.36] claim that spontaneous magnetization occurs in QED 3 . Their analysis is, 
however, based on two-component fermions and assumes a heavy fermion mass which breaks parity 
explicitly. 
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APPENDIX A: CONTOUR INTEGRALS 



In this appendix we provide formulas for contour integrals which appear in Section |l 1 
and Section fV[ They are given for an arbitrary positive constant M and an integer n by 

^< M >4£/*(?rW' (A1) 

T {n) , - ~ It, iz 



^■■=-Jc» dz WTWy' (A2) 

J(°)(M) := - £ dz \nl 2 (z 2 + M 2 ), (A3) 

where the contour C M is a rectangle which connects two trajectories (—R — ifx, +R — iji) and 
(-R, +R), and finally we take the limit R — > oo (see Fig. [[]). 

When = M, poles or branch points of the integrand are put just on the contour C^. 
We therefore replace the integrals in \y\ — M as their Cauchy's principal values. Notice that 
integrals with n > 2 show divergence at = M and become ill-defined. To avoid these 
ambiguities we define them as n-th M derivatives of the integrals I± , namely, 

/ f) (M) = rr ( — I { l\M). ( A4) 

± y ! (n-l)\ \2MdM) ± y ! y ' 

The integrals I±^ are calculated by means of the residue theorem. When > M, the poles 
z = ±iM of the integrand are located inside the contour C^. Thus we obtain the following 
results 

lM(M) = -±e(\n\-M), (A5) 
(M) = — sgn(/i) 9(\fi\ — M), (A6) 



r(i) 



where the step functions take a value 1/2 due to the redefinition as their Cauchy's principal 
value at = M and they match the zero temperature limit of Fermi-Dirac distribution 
function. 

By use of Eq. (|A4p the integrals I± ^ are determined as 

4 2) (M) = - M) - ^ 5(|/z| - M), (A7) 

J - 2)(M) = -^ s e n 0") W - M )' ( A8 ) 

which show delta function singularities at = M as mentioned above. They cannot be 
neglected because they appear as the zero temperature limit of the M derivative of Fermi- 
Dirac distribution function. 

As to the integral 7'°', its integrand does not have any pole but branch points at z — 
±iM. When > M, the branch points are located inside the contour therefore the 
integral is replaced with the one integrated along a minimal contour that rounds the branch 
cut \lmz\ > M on the imaginary axis. Thus we obtain 

I^\M) = 2(|/i| — M) 9(\fx\ — M). (A9) 
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APPENDIX B: ASYMPTOTIC EXPANSION 



In the following we provide the asymptotic expansion for the function I(X) which is 
given by the integral 

I(X) := / dr G(r)e- X ^ T \ (Bl) 
Jo 

We suppose that the function $(r) is analytic at r = 0, as well as the function G(t), and an 
odd function for r which monotonically increases in r > 0. For convenience, we normalize 
$ as 



= 1. (B2) 

T = 



Our aim is to expand I{X) into the asymptotic series for the large X, that is 

OO T 

(* - oo). (B3) 

We notice that there exists the inverse function f(y) := which is also an odd 

function of y and monotonically increases in y > from the definition of $(t). We exchange 
the variable of integration with y = $(r) in Eq. ( |B1| ) and obtain 

I(X) = fdy^ G{f{y)) e -x y {m) 

Because of an exponential factor e~ Xy , the neighborhood of y = dominantly contributes 
to the integral for the large X. It is therefore allowed in Eq. ( B4|) to expand 

F(y):=^G(f(y)) (B5) 
into a power series of y. Thus the asymptotic series of I{X) is obtained as 

W^E^. (X^oc), (B6) 

where denotes the n-th derivative of F(y) at y = 0. The coefficients F^ are determined 
by derivatives of G and /, and the series up to 0(X~ A ) becomes 

r/v , G(°) G (1) G( 2 ) + /( 3 )G(°) G^+3f^G^ 

J(X)~^ + ^ + ^ + ^ + -.., (X-oo). (B7) 

where G^ (J^) denotes the n-th derivative of G (/) at r = (y = 0). 

Now we apply the above expansion to the vacuum polarization functions given in Section 



rV| . The functions Il e , II j_, and II are written in the same form as Eq. ( pl|) , that is, 
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n e (p) = 
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/ dv 
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JO 
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(B8) 



where we define the functions and $(r) by 



:= Z 2 + M 2 , Af„ := \ m 2 + 



1 — f 2 



p 2 , 



$(r) :=r 1 + 



X(/c — i/i) 



cosh r — cosh rf 1 — f 2 



2r sinh r 



cos 2 } , 



(B9) 



where x := (/p) 2 and 9 denotes the polar angle for the Euclidean momentum p^. The 
functions G e , G±, and G Q are read from Eq. ( 4. ICQ as 



GJt) 



Gdr) 



r(cosh tv — v coth r sinh tv ) 

sinh r 
2r(cosh t — cosh rf ) 

sinh 3 t 



G e (p), 



(BIO) 



_ . . 2r(coshr coshrf — 1) 

G (r) = —2 sgn( eJ B), 

smn r 



all of which are analytic at r = 0. 

We notice that $(r) is a monotonically increasing odd function of r and satisfies the 
normalization condition of Eq. (|B2j ). Therefore the vacuum polarization functions are ex- 
panded into the same form as the asymptotic series in Eq. (fBTJ) in the weak field regime 

(Z- 2 «f 2 ); 



n e (p) = 


a/ 2 

7T 


Trff 
JO J 


POO 

/ djfe 

— oo 


n±(p) = 


a/ 2 


Trff 


POO 


7T 


Jo ^ 


— oo 



X(k-i/j) X 2 (k-in) 

Li U I 



a/ 2 f 1 f 00 

n o (p) — — sgn(e£>) dv dk i(k — iji) 
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X(k-ifi) X 2 (k-i/j,) 
Gl 0) 



0(X~ 3 ) 

o(x- 3 ) 



(Bll) 



+ 



X(k-ifi) X 2 (k-i/2) 



Q(X 



-3^ 



Since the only dimensionful parameter is the momentum p after taking the symmetric limit 
J± — > 0, the X _n terms in the above equations are identified with the n-th terms of the 
asymptotic series in Eq. (|4.27[ ), whose coefficients up to 0(X~~ r ) are given by 



pi 1 ("OO ] /■ 1 

= limp (1 - f 2 ) - / — + i cfe — 

If iz 
= lim p 2 sgn( eJ B) / df (1 + f 2 ) -f dz 

J±-*0 JO 7T JC M [2 2 + M 2 ] 



(B12) 



= 0. 
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The other coefficients up to 0(X^ 1 ) become trivially zero due to even or odd property of 
the functions G e , G±, and G Q under the parity. The integrals over z are carried out as the 
contour integrals given in Appendix [A| and are given by 



1 dv 1 —^[l-9(\^\-M v )} 



Hf ] = limp/ 

J±-+o Jo 

H^ 1 ' = — lim p 2 sgn( fieB) / dv 
J±->o Jo 



1 + v 2 
2M„ 



6{\v\-M v ) 



(B13) 



which have appeared in Section |V| and the integrals over parameter v can be calculated by 
using formulas in Appendix |C|. Thus we obtain 



Hf ] = lim \6(m- \fi\)I{p;m) + 9(\fj,\ - m) 6\p 2 - 4(> 2 - m 2 )]l(p; /i) 

J± — > L L J 



H^ 0) = - lim 8ga(jieB) 9(\/j\ - m) e\p 2 - 4(/x 2 - m 2 ) 



4 \p 2 - 2(/i 2 - m 2 )] 



with the function I(p; u) which is defined as 



J(p; u) :-- 



2\u\yp 2 — 4(u 2 — m? 

p2 



+ 1 



4m 
p2 



arctan 



pjp 2 — 4(/i 2 — m 2 ) 



\jp 2 — 4(w 2 — m 2 ) 
2M 



(B14) 



(B15) 



APPENDIX C: FEYNMAN INTEGRALS AT FINITE DENSITIES 



In this paper we have used the following integral formulas in calculations of loop dia- 
grams; 



4(p) 

MP) 



p f 1 dv 



2 Jo M™ 
p f 1 dv 



MP) ■= o 



2 Jo Mr M ~ Mv) 
9(\/i\ - m){l n (m) - 6[p 2 - A(fi 2 - m 2 ) 
p f 1 dv 



S(\fi\ - My 



2 Jo M r v 

= 6{\n\ - m)6[p 2 - 4(/i 2 - m 2 )] /„(//), 
where n is an arbitrary integer and the functions I n and I n are given by 



. , <C(p) - u 2 (l n 3 jtfj 

nW '~ [C(p)YI 2 2 \V 2' 2' C(p), 
i 

:= 



l«l n_ V c(p) - 

with the hypergeometric function 2-^1(0, b, c; z) and 



(CI) 



(C2) 



(C3) 



(C4) 
(C5) 
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C(p) :=m 2 + P - 



For instance the vacuum polarization function II(p) in Eq. ( |6.28 ) is given by 



flip) = a f dv l —^[l-9{\Li\-M v )] 
Jo M v 



8a 
P 3 



{^(p) - m 2 J°(p) - [J_!(p) - m 2 Ji(p)] } 
j#(m — J_i(m) — m 2 /i(m) 



+ 0(|//| -m)#[p 2 -4(/i 2 -m 2 )] [/-i(^) -m 2 /^)]}. 
The functions I_i and Ji are calculated as 



7 _ iM = HVgM^! + gM arctan 



C(p) - u 2 



ii(w) = arctan 



C(p) - it 2 



Finally we obtain 



fl(p) = _{0( m _ |/i|)/(p;m)+ 0(M -m)^[p 2 -4( / u 2 -m 2 )]/(p;/i)}, 
where we define the function J(p; w) as 



I{p;u) := — I-i(u) -m 2 Ii{u) 
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2\u\^Jp 2 — 4(w 2 — m?) 



P 



+ 1 
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FIG. 1. The contour of integration, C M , on the complex z-plane (for > 0). 
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FIG. 2. Possible approaches to the symmetric limit on the (|//|, to) -plane: I. | | < m — > 0, II. 
m < — > 0, and III. = to — > 0. 
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FIG. 3. The Schwinger-Dyson equation. The lines with blob mean the full propagators. 




FIG. 4. The two-particle irreducible (2PI) diagram at the leading order in 1/N expansion. 
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FIG. 5. The b dependence of the function K(x). 
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TABLES 



\k\ 






m d 


iMdl 




a 
2a 


-N/2tt1 2 
-N/AttI 2 





N/Anl 2 
N/2tt1 2 


LO s /l 

u s /2l 





u s /2l 

OJs/l 



TABLE I. Condensates and solutions to the Schwinger-Dyson equation for each k. 
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